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Search for New Physics in the B° H mixing phase at CDF 
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The Collider Detector at Fermilab (CDF) experiment performed the first measurement of the time-evolution of flavor- 
tagged — > J/ipifi decays, which probes mixing-induced CP- violation in the sector. Any sizable deviation from 
zero of the phase (3$^^ , accessible through interference of the b — > ccs quark-level process accompanied or not by 
B® — mixing, would be unambiguous indication of physics beyond the Standard Model. I report CDF results 
obtained in 1.35 fb — 1 , a recent extension to a larger dataset corresponding to 2.8 fb _1 , and future projections. 
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<N . 1. INTRODUCTION 

> 

Many precise results from several years of successful B-factories' running disfavor significant (> 10%) contributions 



from "New Physics" (NP) in tree-dominated 6-meson decays. Agreement with the Standard Model (SM) is also 
found in higher-order processes such as K°-K or B°-B transitions due to second-order weak interactions (mixing) 
that involve virtual massive particles and may receive contributions from NP. A less clear picture is available for 
the B® system. The strength of NP contributions in B° s -B s mixing is constrained by the precise measurement of 
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the oscillation frequency [lj, which disfavors large magnitudes of NP amplitudes. However, knowledge of only the 
frequency leaves the phase of the mixing amplitude unconstrained. Indeed, possible large NP phases are currently 
' not excluded. The mixing phase is accessible through the time-evolution of B® — > J/ip4> decays, which is sensitive 
. to the relative phase between the mixing and the b — > ccs quark-level transition, 0s = Pg M + Pf F - Such 
phase is responsible for CP-violation and is /3f M = arg(— V ts V t * b /V cs V* b ) ss 0.02 in the SM 0]; any sizable deviation 
from this value would be unambiguous evidence of NP If NP contributes a phase (0f F ), this would also enter 
jJ/iji4> _ ^sm _ 20^ p , which is the phase difference between mixing and decay into final states common to B° and 



B° s , and is tiny in the SM: (j) S s M = arg(-M 12 /T 12 ) « 0.004 j^. The phase enters the decay- width difference 

between light and heavy states, AT = T L -T H = 2\T 12 \ cos^j^), which is Ar SM « 2|r 12 | = 0.096 ± 0.036 ps" 1 
in the SM [4| and plays a role in B® — > J/ip4> decays. Since the SM values for J J^"^ and § J J^"^ cannot be resolved 
with the resolution of current experiments, the following approximation is used: w — 2/3f p ss —20s^^, which 

OO , holds in case of sizable NP contributions. 

This measurement of 0^^ is analogous to the determination of the phase (3 — arg(— V c dV* b /VtdV t * b ) in B° — > 
J/ipKg decays, except for a few additional complications: the oscillation frequency is about 35 times higher in 
than in B° mesons, requiring excellent decay-time resolution; the decay of a pseudoscalar meson (B° ) into two vector 
mesons (J/4> and 0) produces two CP-even states (orbital angular momentum L = 0,2), and one CP-odd state 
(L = 1), which should be separated for maximum sensitivity; and the value of the SM expectation for 0s is 
approximately 30 times smaller Q than the known (3 value 



2. SIGNAL SELECTION 

The CDF experiment at the Fermilab Tevatron performed the first measurement of the time-evolution of flavor- 
tagged B® — > J/ip(—> /i + fi~)(f>(-^ K + K~) decays [7]. These were reconstructed in pp collision data corresponding 
to a time-integrated luminosity of 1.35 fb _1 . Events enriched in J ftp decays are selected by a trigger that requires 
the spatial matching between a pair of two-dimensional, oppositely-curved, tracks in the multi-wire drift chamber 
(coverage \rj\ < 1) and their extrapolation outward to track-segments reconstructed in the muon detectors (drift 
chambers and scintillating fibers). In the offline analysis, a kinematic fit to a common space-point is applied between 
the candidate J/ip and another pair of tracks consistent with being kaons originated from a <j) meson decay. The 
measurement of specific energy loss by ionization in the drift chamber (dE / dx) provides 1.5cr separation between 
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Figure 1: Mass distribution with fit projection and sideband regions overlaid (a). Distribution of — 2Aln(L p ) observed in 
pseudo-experiments (solid black line) compared with the nominal \ 2 distribution (solid green line) (b) . The effect of sampling 
the nuisance parameters within 5a of their estimates in data is shown by the dashed colored lines. The vertical lines indicate the 
different values of — 2Aln(L p ) to be used in the ideal x 2 case (green), in the observed distribution (black), and in the observed 
distribution including systematic uncertainties (red), to obtain a 95% C.L. region. The distribution shown corresponds to a 
specific {fiV^^ Ar) point, but similar distributions have been observed across the whole space. 



charged kaons and pions with momenta p > 2 GeV/c. At lower momenta, scintillators bars surrounding the chamber 
measure arrival times of charged particles (time-of-fhght, TOF) with approximately 110 ps resolution, providing 
separation between kaons and pions in excess of 2a. An artificial neural network trained on simulated data (to 
identify signal, S) and B® mass sidebands (for background, B) is used for an unbiased optimization of the selection. 
The quantity S/ \/S + B is maximized using kinematic and particle identification (PID) information. Attempts of 
using the average statistical resolution on f^J^ observed in ensembles of pseudoexperiments as figure of merit were 
inconclusive because of irregularities of the likelihood (see below). Discriminating observables include kaon-likelihood, 
from the combination of dE / dx and TOF information; transverse momenta of the B® and </> mesons; the K + K~ 
mass; and the quality of the vertex fit. The final sample contains approximately 2000 signal events over a comparable 
background (fig. [T] (a)). Seven layers of silicon sensors extending radially up to 22 cm, and the drift chamber that 
provides 96 measurements between 30 and 140 cm, all immersed in the 1.4 T axial magnetic field, provide a mass 
resolution of approximately 10 MeV/c 2 on the B® — > J/ip(f> peak. 



3. FITTING THE TIME EVOLUTION 

The sensitivity to the mixing phase is enhanced if the evolution of CP-even eigenstates, CP-odd eigenstates, and 
their interference is separated. CDF uses the angular distributions of final state particles to statistically determine 
the CP-composition of the signal. The angular distributions are studied in the transversity basis, which allows a 
convenient separation between CP-odd and CP-even terms in the equations of the time-evolution [8j. 

Sensitivity to the phase increases if the evolutions of bottom-strange mesons produced as B® or B s are studied 
independently. The time development of flavor-tagged decays contains terms proportional to sin(2/3/^^), reducing 
the ambiguity with respect to the untagged case (oc | sin(2/3/^^)|). Building on techniques used in the B® mixing 
frequency measurement [l|, the production flavor is inferred using two classes of algorithms. Opposite-side tags 
exploit bb pair production, the dominant source of 6-hadrons at the Tevatron, and estimate the production flavor 
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from the charge of decay products (e, [i, or jet) of the 6-hadron produced from the other 6-quark in the event. Same- 
side tags rely on the charges of associated particles produced in the fragmentation of the 6-quark that hadronizes into 
the candidate B® meson. The tagging power, eD 2 « 4.5%, is the product of an efficiency e, the fraction of candidates 
with a flavor tag, and the square of the dilution D = 1 — 2w, where w is the mistag probability. Multiple tags, if 
any, are combined as independent. The proper time of the decay and its resolution are known on a per-candidate 
basis from the position of the decay vertex, which is determined with an average resolution of approximately 27 /jm 
(90 fs _1 ) in B® — > J/ip4> decays, owing to the first layer of the silicon detector at 1.6 cm radius from the beam. 

Information on B® candidate mass and its uncertainty, angles between final state particles' trajectories (to extract 
the CP-composition), production flavor, and decay length and its resolution are used as observables in a multivariate 
unbinned likelihood fit of the time evolution that accounts for direct decay amplitude, mixing followed by the decay, 
and their interference. Direct CP-violation is expected to be small and is not considered. The fit determines the 
phase f3s ^, the decay-width difference Ar, and 25 other "nuisance" parameters (is). These include the mean B® 
decay-width (r = (r^ + Th)/2), the squared magnitudes of linear polarization amplitudes (|^4o| 2 , I^||| 2 j I^jJ)' * ne 
CP-conserving ("strong") phases ($h = axg(A«Ao), 6± — arg(Aj_A;5)), and others. The acceptance of the detector 
is calculated from a Monte Carlo simulation and found to be consistent with observed angular distributions of 
random combinations of four tracks in data; the angular-mass-lifetime model was validated by measuring lifetime 
and polarization amplitudes in 7800 B° — > J/tpK*° decays, which show angular features similar to the B° s sample: 
ct(B°) = 456 ± 6 (stat.) ± 6 {system, \A \ 2 = 0.569 ± 0.009 (stat.) ± 0.009 (syst.), \A\\ \ 2 = 0.211 ± 0.012 (stat.) ± 
0.006 (syst.), $|| = -2.96±0.08 (stat.) ±0.03 (syst.), and$ ± = 2.97±0.06 (stat.) ±0.01 (syst.). The results, consistent 
and competitive with most recent B-factories' results Q , support the reliability of the model. Additional confidence 
is provided by the precise measurement of lifetime and width-difference in untagged B® — > J /^4> decays [l^] • 



4. STATISTICAL ISSUES 



Tests of the fit on simulated samples show biased, non-Gaussian distributions of estimates and multiple max- 
ima, because the likelihood is invariant under the transformation T = (2/3/^^ — * 7r — 2/?/^^, Ar — > — Ar,$ii — > 
27T — $|i, $l — * 7T — $j_), and the resolution on fii'^ was found to depend crucially on the true values of (3s and Ar. 
CDF quotes therefore a frequentist confidence region in the (/9/ , Ar) plane rather than point-estimates for these 
parameters. Obtaining a correct and meaningful region from a multidimensional likelihood is challenging: one should 
construct the full 27-dimensional region, a difficult task computationally, and project it onto the ((3 J S ^^ ', Ar) plane. 
The choice of the ordering algorithm is critical to prevent the projection from covering most of the (/3 s W ,Ar) space, 
yielding a scarcely informative result. A common approximate method is to replace the likelihood, L(/3s^^, AT, if), 
with the profile likelihood, L p (0^, AT, if). For every point in the (/3s , AT) plane, v are the values of nuisance 
parameters that maximize the likelihood. Then — 2Aln(L p ) is typically used as a \ 2 variable to derive confidence 
regions in the two-dimensional space (pi^^ , AT). However, the simulation shows that in the present case the approx- 
imation fails: the resulting regions contain the true values with lower probability than the nominal confidence level 
(C.L.) because the — 2 A hx(L p ) distribution has longer tails than a \ 2 > an d is not even independent of the true values 
of the nuisance parameters (fig.[T](b)). A full confidence region construction is therefore needed, using simulation of a 
large number of pseudo-experiments to derive the actual distribution of — 2Aln(L p ), with a potential for an excessive 
weakening of the results from systematic uncertainties. However, in a full confidence limit construction, the use of 
— 2Aln(L p ) as ordering function is close to optimal for limiting the impact of systematic uncertainties ll|. With 
this method, CDF is able to rigorously account for the effect of systematic uncertainties just by randomly sampling 
a limited number of points in the space of all nuisance parameters: a specific value (fts^^, AT) is excluded only if it 
can be excluded for any assumed value of the nuisance parameters within 5er of their estimate on data. The result 
is a Ps^^-AT contour that is the truly two-dimensional projection of the full, 27-dimensional confidence region. 
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CDF Run II Preliminary L = 2.8fb 1 
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Figure 2: Confidence region in the (/3 s J/,/ " /> , Ar) plane obtained with 2.8 fb _1 of CDF data (a). The green band is the region 
allowed by any NP contribution not entering |ria|, and assuming 2|Fi2| = 0.096 ± 0.036 ps -1 Q. Fraction of CDF-equivalent 
experiments that would observe a 5cr deviation from the SM as a function of the value of fti^^ in a sample corresponding to 
6 fb (black line) or 8 fb -1 (red line) of integrated luminosity (b). 



5. RESULTS 

The results on 1.35 fb -1 show a 1.5a fluctuation with respect to the SM values Considering Ar as an additional 
nuisance parameter, the 68% C.L. allowed region for the mixing phase is 0.16 < ^ 3 J^ < 1.41, which restricts to 
pJ/M e [0.12,0.68] U [0.89,1.45] assuming no NP contributions in \T 12 \ (i. e., constraining 2|r 12 1 to 0.096 ± 0.036 
ps -1 [4]). This result has been confirmed by the D0 Collaboration [ljj . which observed a consistent fluctuation in an 
analysis where the two-fold symmetry of the likelihood is removed by_ assuming an additional theoretical constraint 
between strong phases of — » J/tpcj) an d B° — > J/ipK*° decays [13J. After removing this assumption, CDF and 
D0 results can be combined yielding a 2.2<r fluctuation with respect to the SM and the following 68% C.L. range: 
ptm e [0.24,0.57] U [0.99, 1.33] Q. 

CDF has reported at this conference a partial extension of the analysis to a larger sample, corresponding to 2.8 
fb -1 . This is approximately equivalent to 2.0 fb -1 effective luminosity, because the calibration of dE / dx and TOF 
was unavailable for the whole sample and PID information is not used in the selection, nor in flavor tagging for the 
second half of the dataset. More than 3200 decays are reconstructed, but approximately 4000 are expected when 
PID will be available in the selection. Figure [2] (a) shows the results. The two regions symmetric with respect 
to the (7r/4, 0) point reflect the symmetry of the likelihood, which cannot determine from data if cos(5j_) < and 
cos(<5^ — #|i ) > (corresponding to the Ar > solution) or viceversa (AT < 0). The fluctuation with respect to the 
SM is confirmed and strengthened, reaching the 1.8<r level. The updated analysis restricts the allowed regions for 
the phase to the range 0.28 < (3^^ < 1.29 at the 68% C.L. 

Although the observed deviations are not yet significant, the pattern of independent results showing consistent 
fluctuations in the same direction is promising in view of the analysis of the full dataset, expected to reach approx- 
imately 6 fb -1 by year 2009, or 8 fb -1 by 2010, if Run II of the Tevatron will be extended. Figure [5] (b) shows 
the probability of a 5tr exclusion of the SM at CDF as a function of the value of fti^^ in these two scenarios and 
assuming Ar = 0.1 ps -1 . This extrapolation, which assumes no external constraints and no improvements in the 
analysis, is conservative: CDF is improving the analysis, with significantly increased tagging power, a 50% additional 
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signal collected by other triggers, and the possibility to resolve the strong-phases ambiguity using data [14(; tight 
external constraints (e. g., on the B® lifetime) can be applied, and CDF and D0 results will be combined for maxi- 
mum Tevatron sensitivity. As happened in the past, deviations from expectations in measurements of lower-energy 
processes may indicate NP prior to direct discovery of new resonances, as those expected in the forthcoming run of 



the Large Hadron Collider 15 1. 
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